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Abstract 

We determine the allowed parameter space and the particle spectra of the minimal SUSY 
standard model (MSSM) and nonminimal SUSY standard model (NMSSM) imposing 
correct electroweak gauge-symmetry breaking and recent experimental constraints. The 
parameters of the models are evolved with the SUSY RGEs assuming universality at the 
GUT scale. Applying the new unbounded from below (UFB) constraints we can exclude 
the LSP singlinos and light scalar and pseudoscalar Higgs singlets of the NMSSM. This 
exclusion removes the experimental possibility to distinguish between the MSSM and 
NMSSM via the recently proposed search for an additional cascade produced in the decay 
of the bino into the LSP singlino. Furthermore the effects of the dark matter condition 
for the MSSM and NMSSM are investigated and the differences concerning the parameter 
space, the SUSY particle and Higgs sector are discussed. 



1 Introduction 



If the standard model (SM) is embedded in a grand unified theory (GUT), finetuning is 
needed in order to have a Higgs mass of the order of the electroweak scale and not the 
GUT scale. The quadratic divergences in the radiative corrections to the Higgs mass are 
the reason for this naturalness problem in the SM. SUSY models have the advantage of 
being free of quadratic divergences due to cancellations in Feynman graphs with particles 
and their supersymmetric partners. The simplest SUSY extension of the SM is the mini- 
mal SUSY standard model (MSSM) [[1], || with two Higgs doublets Hi and H 2 . However 
the MSSM has also a naturalness problem why the /i-parameter of the term \i Hi H 2 
in the superpotential should be much smaller than the Planck mass. The nonminimal 
SUSY standard model (NMSSM) is a minimal extension of the MSSM. In this model 
a Higgs singlet N is added to the Higgs sector of the MSSM. Moreover the /i-parameter 
of the MSSM can be dynamically generated via \i = Xx with the term A Hi H 2 N in 
the superpotential of the NMSSM, if the Higgs singlet acquires a vaccum expectation 
value (VEV) (N) = x. SUSY models with Higgs singlets can be derived from superstring 
inspired Eq or SU(5) x U(l) GUT models and they offer the possibility of spontaneous 
breaking of the CP symmetry. The discrete Z 3 symmetry of the NMSSM causes a domain 
wall problem. In the early universe the Z3 symmetry is spontaneously broken during the 
electoweak phase transition. Domains of different degenerate vacua develop which are 
separated by domain walls. Since these domain walls would dominate the energy density 
of the universe and disturb primordial nucleosynthesis, they have to disappear before the 
onset of nucleosynthesis. Possible solutions to the domain wall problem and the hierarchy 
problem reintroduce the terms (i Hi H 2 and // N 2 in the superpotential of the NMSSM 
and possess a gauged R-symmetry or a target space duality symmetry at the Planck scale 
||. In these more general models the NMSSM is approximated for very small values of 
the parameters \i and //. 

The addition of a Higgs singlet superfield leads to two extra Higgses and one extra neu- 
tralino in the NMSSM compared to the particle spectrum of the MSSM. Since the Higgs 
singlets and the singlino N can mix with the other Higgses and neutralinos, the phe- 
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nomenology of the NMSSM may be modified compared to the MSSM. The phenomenol- 
ogy of the NMSSM has been investigated in several papers - In SUSY models 
with R-parity conservation the lightest SUSY particle (LSP) cannot decay and is therefore 
stable. Especially a neutral LSP like the lightest neutralino is a good candidate for cold 
dark matter. In the MSSM the dark matter neutralino is mostly a bino. In the NMSSM 
there exists the alternative of a dark matter neutralino with a larger singlino portion. 



The cosmology of these dark matter singlinos has been discussed in |16| at a given fixed 
low energy scale. The experimentally and cosmologically allowed parameter space of the 
NMSSM and the dark matter neutralinos have been investigated in |T7| using SUSY RG 
evolutions. 

In the present paper we compare the allowed parameter space and the particle spectra of 
the MSSM and NMSSM. Assuming universality at the GUT scale the parameters of the 
models are evolved to the electroweak scale with the SUSY RGEs (given in [|18| and Ap- 
pendix A2 for the NMSSM). The Higgs potential is minimized and additional theoretical 
constraints and recent experimental bounds are imposed. Besides charge and color break- 
ing (CCB) constraints we include here new unbounded from below (UFB) constraints 
Furthermore the consequences of the dark matter condition for the MSSM and NMSSM 
are investigated. In the present paper we give a more detailed analysis of the two models 
and the dark matter condition than in our previous paper ||20|| . 

In section 2 we introduce the NMSSM. The additional theoretical and experimental con- 
straints are discussed in section 3. In section 4 we describe the dark matter constraints. 
The procedure of finding solutions of the NMSSM is explained in section 5. In section 
6 the parameter space and the particle spectra of the MSSM and NMSSM are investi- 
gated without and with the dark matter condition. In Appendix Al we give the complete 
formulas for the neutralino annihilation cross section in the NMSSM which are needed 



to calculate the relic neutralino density. In Appendix A2 the SUSY RGEs |1| for the 
NMSSM are listed. 
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2 The NMSSM 



The Higgs sector of the MSSM consists of two Higgs doublets H\ and H 2 . In the NMSSM 
H an additional Higgs singlet N is introduced. The superpotential of the NMSSM con- 
tains only Yukawa terms with dimensionless couplings 

W = h d Hl eQ D - h u Hl eQU + X Hf e H 2 N - - kN 3 . (1) 

3 

Q is a squark doublet. U and D are up-type and down- type squark singlets, e is the 
antisymmetric tensor with e\ 2 = 1. The SUSY breaking is parametrized by different 
types of SUSY soft breaking terms in the lagrangian: 

-h d A d H^eQD + h u A u eQU 
+ XA X H? e H 2 N +\kA k N 3 

+ 1 M 3 A 3 A 3 + ~ M 2 A* A" + i M 1 A a \i + h.c. (2) 
In SUSY models the tree-level scalar potential is given by the formula: 

V tree = \ E(D a ) 2 + J2 F i F i + ^soft (3) 

D a = g a A\T^A h (4) 

«=-m <«> 

where A{ are scalar fields and T°j are the gauge group generators. The part of the scalar 
potential with only Higgs fields corresponds to the Higgs potential. The neutral part of 
the tree-level Higgs potential is given by: 

V*2» = m 2 Hi l^l 2 + m 2 H2 \H° 2 \ 2 + m 2 N \N\ 2 + A 2 (l^ 2 + |tf 2 °| 2 ) |iV| 2 + A 2 |^ | 2 |tf 2 °| 2 
+ k 2 |iV| 4 — Ak N* 2 — X A\ Hi H°N-^kA k N 3 

+ l(9 2 y + 9l)(\H° 1 \ 2 -\H°\ 2 ) 2 + h.c. (6) 

In order to calculate the Higgs VEVs reliably we use the effective 1 loop Higgs potential 
consisting of the tree-level Higgs potential and the radiative corrections 

y Higgs _ y Higgs . y Higgs /y\ 
Hoop tree ' v rad " V / 



The radiative corrections |J to the effective Higgs potential arising from all particles and 
antiparticles with field-dependent mass m,j, spin Si and color degrees of freedom C, are 
given by 



Q is the renormalisation scale which we take to be equal to the electroweak scale M wea k — 
100 GeV. Here we only consider the important contributions of the top quark and stops 
ti t 2- The electroweak gauge-symmetry SU(2)j x U(l)y is spontaneously broken to the 
electromagnetic gauge-symmetry U(l) em by the Higgs VEVs (Hf) = Vi with i — 1,2 and 
(N) = x. The local minimum of the effective 1 loop Higgs potential is determined by the 
three minimum conditions for the VEVs: 

1 2 m 2 Hl + S 1 - (m 2 H2 + E 2 ) tan 2 /? 2 

2™* = ta^^l AX (9) 

-• (9R\ - 2\x(A x + kx) 

Sm ™ ~ + S 1 + m? H2 + T 2 + 2 A 2 x 2 + A 2 v 2 ( } 

(m^ + S 3 ) x 2 - k A k x 3 + 2 k 2 x 4 + A 2 x 2 v 2 - ~ (A A + 2 fc x) A x sin(2/?) = (11) 

where u = ^/t> 2 + v 2 = 174 GeV, tan/3 = v 2 /v 1 and S ? = dV r Jf 9s / dv 2 . 

The neutralinos are Majorana particles and mixed states of the gauginos and Higgsinos. 

The mass and composition of the neutralinos is defined by the following part of the 

lagrangian 

£ = -^ T M^ + h.c. (12) 

V T = {-iXx, K^K^n). (13) 

In this basis the symmetric mass matrix M of the neutralinos has the form: 

/ Mi — m z sin 6 W cos (3 m z sin 6 W sin (3 

M 2 m z cos 6 W cos f3 —m z cos 6 W sin (3 

— m z sin 6w cos (3 m z cos ^ cos (3 Ax A i> 2 

m z sin 6^ sin /3 — cos sin (3 Ax A v\ 

Aw 2 A«i -2A;x 



4 



With fj, = A x the first 4x4 submatrix recovers the mass matrix of the MSSM 0. The 
mass of the neutralinos is here obtained by diagonalising the mass matrix M with the 
orthogonal matrix N. (Then some mass eigenvalues may be negative pl| .) 

C = -\m l $ix Q i (14) 

,,o ( X° 



X- = ( ^ J with x- = NyVj and M diag = N M N T . (15) 

The neutralinos Xi {i = 1 — 5) are ordered with increasing mass \uii\, thus Xi is the 
LSP neutralino. The matrix elements = 1 — 5) describe the composition of 

the neutralino Xi m the basis \&j (Eq.(|l"3D). For example the bino portion of the LSP 
neutralino is given by Nf x and the singlino portion of the LSP neutralino by Nf 5 . The 
diagonalising of the mass matrix M leads to a polynomial equation of degree 5, which 
cannot be solved analytically in general. Therefore we calculate the eigenvalues and 
eigenvectors of the symmetric mass matrix M numerically by using a subroutine in the 
NAG Fortran Library. 

The possibility of gauge coupling unification is an important motivation for SUSY GUTs. 



In contrast to non-supersymmetric GUTs [22], the gauge couplings in SUSY GUTs are 



allowed to unify g a (M x ) = #5 — 0.72 at a scale M x — 1.6 x 10 16 GeV. Within a minimal 
N = 1 supergravity framework universality naturally arises for the SUSY soft breaking 
parameters: 

mi(Mx) = m 

M a (M x ) = m 1/2 (16) 
MM X ) = A (A X (M X ) = -A ). 

The Yukawa couplings at the GUT scale take the values X(M X ) = Ao, k(M x ) = ko 

and h t (M x ) = h t0 . With the SUSY RGEs (@, Appendix A2) the SUSY soft breaking 

parameters and the Yukawa couplings are evolved from the GUT scale to the electroweak 

scale M wea k- At M wea k we minimize the effective 1 loop Higgs potential V^^ s (Eq. ([<])). 

3 Theoretical and experimental constraints 

The solutions of the SUSY RGEs (10, Appendix A2) have to fulfill several theoretical 
and experimental constraints. The minimum of the effective 1 loop Higgs potential V n ^ s 



(Eq.(]7|)) should be not only a local minimum but also a stable minimum. A local minimum 
is stable against tunnelling into a lower minimum, if its lifetime is large enough (of the 
order of the age of the universe). A global minimum is certainly stable, but it may be 
difficult to prove, that it is really a global minimum. Here we only investigate, whether 
the physical minimum is lower than the unphysical minima of V^^ s (vi, v%, x) with at 
least one vanishing VEV 0. 

The conservation of charge and color would be broken by slepton and squark VEVs. To 
avoid a scalar potential with such a minimum, the trilinear couplings have to fulfill the 
following conditions 



A\. < 3(m 2 H2 + m\, + rrifj.) at scale Q ~ A Ui /h Ui (17) 
A% < 3(rr4 + m%. + m 2 Di ) at scale Q ~ A d Jh di (18) 
Al < 3(m 2 Hi + ml. + m|.) at scale Q ~ AJh ei . (19) 

In |19| improved charge and color breaking (CCB) constraints and unbounded from below 
(UFB) constraints are analysed for the MSSM. The by far strongest of these new con- 
straints is the UFB-3 constraint (Eq.(33) in the first publication of Ref.[[19||), in contrast 
to the remaining ones. To apply the UFB-3 constraint to the NMSSM we have to modify 
the formulas in ||19|| : In the NMSSM a dangerous field direction of the tree- level scalar 
potential Vt ree (Eq.(§)) occurs for VEVs of the neutral Higgs H®, the electron sneutrino 
v e denoted by L x = V\ and the staus f£ and denoted by L 3 and E 3 . The VEVs of 
the staus are introduced in order to cancel the F-term of if°: 



dW 2 



\xH° + h T Lz E 3 \ 2 = 0. (20) 



d Hi 

For the VEVs of the staus we choose the relation e = = E 3 . In analogy to the 
MSSM we have replaced in Eq . (|20D // by A x where the Higgs singlet VEV (N) = x is 
determined by the physical minimum of the NMSSM Higgs potential. Since we consider 
here the scalar potential and not only the Higgs potential, the Higgs singlet VEV x may 
be optimized further. But here we don't investigate this more general case. The NMSSM 
tree-level scalar potential along the described field direction is given by: 

t / 2 1 tr0|2 , / 2 1 2 \ | |2 , 2 i r 1 2 

Vtree = m H2 I #2 I + ( m L 3 + m E 3 ) M + m L 1 \ L l\ 



+ \(/v+d) (l^ 2 T + N 2 -l^i 2 ) 2 

2 

+ m 2 N x 2 k A k x 3 + k 2 x 4 . (21) 

3 

Now we can determine the dangerous field direction of V tree (Eq.(|21|)), which is given by: 

I^| 2 = ^ + I^ 2 T + Ie| 2 , (22) 

y y ' U2 



\\x\ 



\m\ (23) 



provided that 



W ' 1 > U Thf + g[T? 2 - 2K- (24) 



For smaller values of {H^l we obtain \L±\ — 0. In the field direction of Eq.(|2"2"D and ( |2"5| ) 
the tree-level scalar potential of the NMSSM reads: 

V UF B-3 = {m 2 H2 + m 2 Li ) | H% \ 2 + ^ {m\ 3 + m% + m 2 Li ) \ \ 



9 2 + 9i " 3 
For \Li\ = the tree-level scalar potential is given by: 



h T 

Ll + m 2 N x 2 --kA k x 3 + k 2 x 4 . (25) 



2 mi Q „ 2 



2 I 1 1 1 2 l^^l / 2 , 2 \ I rjOl 



V UFB - 3 = m 2 H2 \H»\ 2 + K 3 + m^) |if 2 u 



+ ^ + S 2 2 ) (l# 2 T + ^l# 2 °l) +<^ 2 -^A fc x 3 + fc 2 x 4 . (26) 

Since the tree- level scalar potential Vufbs should be larger than the physical minimum 
of the tree-level Higgs potential v££ 9a (Eq.(|)), the UFB-3 constraint can be written as 

V UF bs{Q = Q) > ( V™° ga (Q = M s ) ) Mm . (27) 

The tree-level Higgs potential V t ^ 9S (Eq.(^)) is minimized at the scale Q = Ms, where the 
radiative corrections V^ ggs (Eq.([8D) are very small. Ms turns out to be a certain average 
of SUSY masses. For the same reason the tree-level scalar potential Vufbs is evaluated 
at the scale Q = Q ~ Max(g 2 |e|,g 2 \H%\,h t li?^,*^ \Li\,M s ). Since Vufbs(Q — Q) 
depends on {H^l the relation fl2"7|) should be fulfilled for all [H^l values below Mj. 



Because we assume h t 3> hb,h T , we restrict |tan/3 | to be smaller than 20 
Higgs singlet VEV x we choose the range \x\ < 80 TeV pi 



For the 



13 1 . In the present paper 



we choose the parameters m and |m!/ 2 | to be smaller than 1500 GeV. This restriction is 



motivated by requiring the absence of too much fine-tuning | 13fl , which disfavors solutions 
with enormous sensitivity of mz to the GUT scale parameters h t o or Ao and k$. The 
dark matter constraints further restrict the allowed parameter region of mo and |mi/2| as 
shown in Fig.|l|. 

For the top quark pole mass we take the range 169 GeV < rrit < 181 GeV corresponding 
to the recent world average measurement |23||. Since charged or colored massive stable 



particles are cosmologically disfavored LSPs, we assume the lightest neutralino to be the 
LSP and impose the constraints m^o < m fl and m^o < m^. The SUSY particles in 
the NMSSM have to fulfill the following experimental conditions: > 41.8 GeV 



m, 



> 70GeV f2§, rriq > mGeV(iim~ g < 300GeU) fZ§ and m- q > 70Ge1/(for all 



m k > 63 GeV p5f, m~ g > 173 GeV |7|, m~+ > 85.4 GeV ||, m H+ > 52 GeV 



^Si > lAGeV |^0[ or the coupling of the lightest scalar Higgs Si to the Z boson is reduced 
compared to the SM BT| assuming visible decays with branching ratios like the SM Higgs 
§,and 

YnZ^x^) <30MeV H 



T(Z^x° 1 x° 1 )<7MeV H 
BR{Z - x?X?) < 10- 5 , (1,1) II. 



4 Dark matter constraints 



Inflationary cosmological models are proposed to solve the flatness, horizon and magnetic 
monopole problem. In these models the universe is flat with a density equal to the 
critical density p cr a = 1.054 x 10 -5 /ig GeV/ cm 3 implying a cosmological density parameter 
of Q = p/pcrit = 1- The Hubble constant Hq = 100 ho km s -1 Mpc~ l is restricted 
by observations to the range 0.4 < h < 1 of the scaled Hubble constant h . To be 
consistent with the observed abundances of the isotopes D, 3 He, 4 He and 7 Li the big bang 
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nucleosynthesis model restricts the baryon density to flbaryons < 0.1. The discrepancy 
between the cosmo logical density Q = pj p crit = 1 and the baryon density requires a 
huge amount of dark matter in the universe. One of the favored theories to explain the 
structure formation of the universe is the cold + hot dark matter model (CHDM) f33| . 
In this model the dark matter consists of hot dark matter (massiv neutrinos) and cold 
dark matter (neutralinos) with Qhot — 0.3 and Q co id — 0.65. With the allowed range of 
the scaled Hubble constant ho the condition for the relic density of the neutralinos in the 
CHDM model reads 

0.1 < tt x h 2 < 0.65. (28) 

The relic neutralino density can be calculated by solving the Boltzmann equation for the 
early universe. The Boltzmann equation describes the particle densities in the cosmic 
plasma governed by the annihilation cross section. The annihilation cross section of the 
neutralinos can be expanded in powers of the relative velocity v of the non-relativistic 
neutralinos: o ann v = a + bv 2 . The thermally averaged annihilation cross section p4] - [3H] 



is then given by < cr ann v > = x 1 ' 5 / (2 y/lx) f °° dv v 2 exp{— x v 2 / 4} a ann v = a + 6b / x 
with x = \m^o\ I T. During the expansion of the universe its temperature falls and 
reaches the freeze-out temperature of the neutralinos Tf r , at which the neutralinos drop 
out of the thermal equilibrium with the cosmic plasma. The freeze-out temperature of 



the neutralinos Tf r follows from the equation |34| - [37| 



x fr = In ( 0.0764 \m^o | M P < GannV > c (2 + c) J , (29) 

V ' \J 9* x fr J 

where Mp = 1.221 x lO^Ge^/ is the Planck mass, g* m 81 is the effective number of degrees 
of freedom at Tf r and c = 1/2. With J(xf r ) = dx / x 2 < a ann v >= a/ Xf r + 3b / x 2 r 



34| - ]36|] the relic neutralino density |34| - |37| is given by 

1.07 x 10 9 1 



Qxk2 °~ ^g-M P J{x fr )GeV (30) 



The expansion of the annihilation cross section a ann v = a+bv is not a good approximation 



in the vicinity of poles or thresholds [pq] . In the present work the case of coannihilation 
is not important [[17], [39[, since the LSP neutralino in the MSSM and NMSSM is mostly 
bino-like or a mixed state with small mass (see Fig.f|). Furthermore the LSP singlinos 



H [K| [L3| of the NMSSM are now excluded by the UFB-3 constraint (Eq.fl27|)) as discussed 
in section 6. 

Compared to the MSSM the annihilation cross section of the neutralinos changes in the 
NMSSM because of more particles in this model, which can mix, and modified vertices. In 
our numerical analysis we consider all possible decay channels. The relevant formulas for 
the NMSSM are given in the Appendix Al. They are obtained by substituting the MSSM 
couplings in ]34j by the corresponding NMSSM couplings and by slightly modifying the 
partial wave amplitudes in 



5 Procedure of finding solutions 

The SUSY soft breaking parameters rrii and A, and the Yukawa couplings A, k and h t 
are evolved with the SUSY RGEs ([18[], Appendix A2) from the GUT scale Mx to the 
electroweak scale M wea k- Since we assume universality at the GUT scale Mx (Eq.(|l6|)) 
we have only 6 GUT scale parameters m , m^, A , A , k and h t0 . Together with the low 
energy parameters tan j3 and x there are 8 parameters in the NMSSM with only 5 of them 
being independent. In our procedure we take A , k , h t o, tan j3 and x as independent input 
parameters. The remaining 3 parameters (mo, m^, Aq) are then determined by the three 
minimum conditions (p|-|TT|) and are numerically calculated by the method described in 
|7], 20| : The three minimum conditions (]9|-|iT|) depend only on the low energy parameters 



rrii, Ai, A, k and h t generically denoted by Pk(M wea k): 

F i ( Pk ) = (i = l-3). (31) 



d Vii 00 p 



dvi 

Since Ao, k and h t o are (randomly choosen) fixed input parameters, all low energy pa- 
rameters Pk{M wea k) are functions fk (which can be constant) of the GUT scale parameters 
m o, m i/2 and A only: 

Pk(M weak ) = fk{m ,mi/ 2 ,Ao). (32) 



To obtain the functions fk we numerically solve the SUSY RGEs (JTB|, Appendix A2) 
dpk /dt = hk(t,pi, ...,p n ) {k = 1 — n). A solution to a system of first-order differential 
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equations can be found with the method of Runge, Kutta and Merson. A subroutine 
for this numerical method is available in the CERN Program Library. An additional 
subroutine has to be written which calculates the derivatives dpk/dt (k = 1— n) depending 
on the current values of t and pi, p n . In the three minimum conditions (j3lD we then 
substitute the low energy parameters Pk(M wea k) by the functions fk and obtain the generic 
equations 

i ? i(/fc(%mi/2,4)) = Gi(m ,m 1/ 2,4o) = fl (i = 1 - 3). (33) 
The SUSY breaking parameters mo, mi/2 and Aq are then determined by the zero points of 



the generic equations Gj(mo, mj/2, Aq) = 0(i = 1—3). In order to find a solution of Eq.(33) 
we first minimize ^i=i C«( m o, m i/2, ^4o) 2 and then we check, whether the minimum is 
also a zero point. The minimization is performed by using the simplex method of Nelder 
and Mead. A corresponding minimization routine SIMPLEX is available in the MINUIT 
package of the CERN Program Library. The numerical procedure needs starting values 
for the parameters mo, mi/2 and Aq, which are randomly generated using a logarithmic 
measure within their ranges. 

To obtain a representative sample of solutions of the NMSSM, we randomly generate 
~ 4.3 x 10 8 points in the 5 dimensional parameter space of the input parameters Ao, ho, 
h tQ , tan/3 and x using a logarithmic measure. The SUSY breaking parameters m , mi/2 
and Aq are then numerically calculated with the three minimum equations (|33|). Having 
found a solution of the SUSY RGEs (HI5|, Appendix A2), whose effective 1 loop Higgs 
potential V^^ s (Eq.(^)) has a local minimum for the input parameters tan/3 and x, 
we impose the additional theoretical and experimental constraints described in section 3. 
The very strong UFB-3 constraint (Eq. (pTD) excludes about 80% of the solutions which 
fulfill all the other constraints of section 3. For every solution we calculate the masses 
of all SUSY particles 0, ||| and Higgses || ||, [12|, [13|]. Altogether we obtain about 6600 



solutions which fulfill all requirements. The ranges of the parameters, for which we find 
solutions, are given by m G [ 60 GeV , 1500 GeV }, \m 1/2 \ G [ 120 GeV , 1500 GeV ], 
\A \ G [180GeV, 6000 GeV], A G [2 x 10~ 3 , 0.5], k G [lO^ 3 , 0.6], h t0 G [0.4, 0.9], 
| tan 0| G [ 2.4 , 20 ] and \x\ G [ 800 GeV , 80000 GeV }. Applying the dark matter 
constraint to the 6600 solutions, about 1500 solutions remain which are cosmologically 
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acceptable. 



6 Parameter space and particle spectra 

In the following we investigate the parameter space and the particle spectra of the MSSM 
and NMSSM. The consequences of the dark matter condition for the two models are 
discussed. We show our results in the form of scatter-plots. Each point represents a 
solution of the SUSY RGEs (given in Appendix A2 for the NMSSM) with additional 
theoretical and experimental constraints imposed as described for the NMSSM in section 
3. In the upper left picture of Fig.[I] the MSSM solutions are shown in the (|mx/2 1 3 m o) 
plane. The parameter region with small mo and large |mi/2| is excluded by the UFB-3 
constraint |Tj|. Furthermore in the region |mi/2| ^> thq the stau or stop would be the LSP. 
But a charged or colored LSP like the stau or stop is cosmologically disfavored [f|lj. In 
the upper right picture the cosmologically allowed solutions are shown. The dark matter 
condition leads to an upper bound of 600 GeV for |mi/2| and leaves mo unconstrained. 
Values of m larger than about 300 GeV are cosmologically allowed due to Z-pole or Higgs- 
pole enhanced neutralino pair annihilation. Comparable scatter-plots for the constrained 



MSSM are discussed in In that paper upper bounds of 9A0GeV for |m.i/ 2 | and 500GeV 
for mo are deduced for m t = 170GeV. The lower left picture shows the NMSSM solutions. 
The effect of the UFB-3 constraint (Eq.fl27|)) in the NMSSM is comparable to the MSSM. 
In the NMSSM this constraint also excludes the parameter region with small mo and large 
|mi/ 2 |. The reason, why the case m 3> \mu2\ is suppressed in the NMSSM, is discussed in 
0. It follows from the condition, that the minimum of the scalar potential does not break 
the conservation of charge and color, and from the condition that the physical minimum 
is lower than the symmetric vacuum, 9 m\ < A\. In the lower right picture the dark 
matter condition is imposed leading to an upper bound of 650 GeV for jm^]. In contrast 
to the MSSM the upper bound for m is about 300 GeV in the NMSSM. These upper 
bounds follow from the cosmological upper bounds of 2Q0GeV for the LSP neutralino and 



450 GeV for the lighter selectron as discussed in |20] and below. This can be understood 



from the fact pi], that the LSP neutralino in the NMSSM is mostly bino-like (see Fig.P 
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and the annihilation cross section is dominated by sfermion exchange. 



For the MSSM we show in the upper left picture of Fig.0 the solution points in the 



plane (m , \Aq\). The upper bound for \A \ follows from the UFB-3 constraint [19|. The 



upper right picture shows the cosmologically restricted solutions. The concentration of 
points corresponds to mo values below about 300 GeV in the upper right picture of Fig.[l|. 
Larger values of mo are cosmologically allowed because of Z-pole or Higgs-pole enhanced 
neutralino pair annihilation. The lower left picture shows the allowed values of |^4o| in the 
NMSSM. Comparable with the MSSM the UFB-3 constraint (Eq.(|7D) leads to an upper 
bound for \A \ in the NMSSM. In contrast to the MSSM there is also a lower bound for 
\A \ in the NMSSM. The NMSSM solutions have to fulfill the relation 9 m 2 < A 2 |§, 
which follows from the condition that the physical minimum has to be lower than the 
symmetric vacuum. The dark matter condition strongly restricts the allowed range of 
\Aq\ as shown in the lower right picture. Larger values of tuq are cosmologically excluded, 
since m ^> \my 2 \ is disfavored in the NMSSM as already discussed. 

In the upper left picture of Fig|3] we show the MSSM solution points in the plane (/z, |M 2 |). 
The parameters /i and M 2 are taken at the electroweak scale M wea k. The solutions are 
lying below the diagonals \M 2 \ = The reason is the first minimum condition which 
favors larger values of for | tan/3| — > 1. The first minimum condition of the MSSM has 
the same form as the first minimum condition in the NMSSM (Eq.(£j|)) if /i is substituted 
for A x. This property of the solutions leads to bino-like LSP neutralinos in the MSSM 
(see Fig.|]). In the upper right picture the dark matter condition is imposed. The upper 
bound of 500 GeV for \M 2 \ corresponds to the upper bound of 600 GeV for |mi/2| in the 
upper right picture of Fig.[TJ. The NMSSM solutions in the lower left picture are similar 
to the MSSM solutions, but they are more restricted. This restriction follows from the 
third minimum condition (Eq. (|TTl) ) in the NMSSM which would be an extra constraint 
for the MSSM. The lower right picture shows the cosmologically allowed solutions. The 
upper bound of 550 GeV for \M 2 \ follows from the upper bound of 650 GeV for |m!/ 2 | in 
the lower right picture of Fig.|l]. 

In the upper left picture of Fig.f| the bino portion of the LSP neutralino | < B \ X\ > | 2 
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versus the LSP neutralino mass is shown for the MSSM. The LSP neutralinos are mostly 
binos | < B \ Xi > | 2 ~ 1, which follows from the first minimum condition in the MSSM 
corresponding to Eq.(§) as already discussed. The cosmologically allowed LSP neutralinos 
are shown in the upper right picture. The dark matter condition sets an upper bound 
of 240 GeV for the LSP neutralino mass corresponding to the upper bound of 600 GeV 
for |mi/2| in the upper right picture of Fig.[l]. The lower bound for the mass of the LSP 
neutralinos is about AOGeV in the MSSM. For the NMSSM the lower left picture shows the 
bino portion of the LSP neutralinos. Similarly to the MSSM most of the LSP neutralinos 
are binos in the NMSSM. If the UFB-3 constraint (Eq. (B7T)) is not applied, there exist 
a large fraction of decoupled pure LSP singlinos with | < B \ \\ > \ 2 = and a much 



smaller fraction of mixed LSP singlinos. In JT0[ the phenomenology of the LSP singlinos 



is investigated. The LSP singlinos are now excluded by the UFB-3 constraint (Eq.(|2"7p), 
which excludes the parameter region with small mo and large |mi/2|. This is just the 



appropriate parameter region for LSP singlinos |9|, [10], [13). The lower right picture shows 
the cosmologically allowed LSP neutralinos in the NMSSM. The upper bound of 260 GeV 
for the LSP neutralino mass corresponds to the upper bound of 650 GeV for |mi/2| in the 
lower right picture of Fig.|l[ The dark matter condition improves the lower bound for the 



LSP neutralinos in the NMSSM to 50 GeV as discussed in [g. In the NMSSM all LSP 
neutralinos with mass m^o « m^/2 are cosmologically excluded, since their coupling to 
the Z boson is larger due to a larger Higgsino portion. 

In the upper left picture of Fig|] the mass of the lighter chargino versus the mass of the 
lighter selectron is shown for the MSSM. The upper right picture shows the cosmologically 
allowed lighter chargino masses. The dark matter condition sets un upper bound of 
450 GeV for the lighter chargino and disfavors somewhat selectron masses larger than 
about 1 TeV. Selectrons with 300 GeV < mg fl are cosmologically allowed due to Z-pole 
or Higgs-pole enhanced neutralino pair annihilation [|^]. For the NMSSM the lower left 



picture shows the mass of the lighter chargino versus the mass of the lighter selectron. 
Compared to the MSSM there are not so many solutions in the NMSSM with small 
chargino mass and large selectron mass following from the suppression of the domain 
m \mi/2\ in the NMSSM as discussed in 0. In the lower right picture the dark matter 
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condition is imposed leading to an upper bound of 500 GeV for the lighter chargino and a 
cosmologically allowed mass range of 100 GeV < mg R < 450 GeV for the lighter selectron 
in the NMSSM. The cosmologically allowed solutions in the MSSM and NMSSM show 
some differences. Only in the MSSM there are cosmologically allowed heavier selectrons. 
In contrast to the MSSM, the dark matter condition leads to an slightly improved lower 
bound of 90 GeV for the lighter chargino in the NMSSM. The reason is the enhanced 
lower bound of 50 GeV for the cosmologically allowed LSP neutralinos in the NMSSM 
as discussed and shown in Fig.§ and the proportionality between the mass of the lighter 
chargino and the mass of the bino-like LSP neutralino. 

In the upper left picture of Fig|] the mass of the lighter u-squark versus the gluino mass is 
shown for the MSSM. The upper boundary line for the lighter u-squark mass results from 
the restriction of m . As mentioned in section 3 we have chosen m to be smaller than 
1500 GeV in order to avoid too much finetuning. The lower bound for the gluino and the 
lighter u-squark is 250 GeV. In the upper right picture the cosmologically allowed lighter 
u-squark masses are shown. In many cases the cosmologically allowed lighter u-squark 
mass is proportional to the gluino mass. The deviations from proportionality with large 
u-squark masses follow from the Z-poles or Higgs-poles shown in the upper right picture 
of Fig.l. The dark matter condition leads to an upper bound of 1600 GeV for the gluino 
following from the cosmological upper bound of 240 GeV for the bino-like LSP neutralino. 
In the lower left picture the mass of the lighter u-squark is shown for the NMSSM. In 
contrast to the MSSM [53], the squark masses in the NMSSM are nearly proportional to 



the gluino mass as discussed in p0"| . The reason is, that in the NMSSM solutions with 
mo ^> 1 77ii/2 1 are somewhat disfavored. In contrast to the MSSM the lower bound for the 
gluino is 330 GeV in the NMSSM, which follows from the lower experimental bound of 
85.4 GeV |28[] for the chargino and from the lower bound of about 2.4 for | tan/3| in the 
NMSSM (see Fig.||) as discussed in ||. Excluding small values of m the UFB-3 constraint 
(Eq. (|27|)) further improves the lower bound for the lighter u-squark to 370 GeV. In the 
lower right picture the dark matter condition is imposed. In contrast to the MSSM, the 
dark matter condition improves the lower bound for the gluino to 370 GeV in the NMSSM 



as discussed in 20 . The reason is the enhanced lower bound of 50 GeV and the larger 
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Higgsino portion of the lightest cosmologically allowed LSP neutralinos in the NMSSM. 
With the dark matter condition the lower bound for the lighter u-squark slightly improves 
to 390 GeV. The cosmological upper bound of 1800 GeV for the gluino in the NMSSM is 
similar to the corresponding upper bound in the MSSM. The upper bound of 1600 GeV 
for the lighter u-squark in the NMSSM corresponds to the upper bound for the lighter 
u-squark in the MSSM at large gluino masses. 

In the upper left picture of Fig]7] the mass of the lighter stop versus the gluino mass is 
shown for the MSSM. The stop masses are restricted from above, since the mo values 
have been choosen to be smaller than 1500 GeV as discussed in section 3. In the upper 
right picture the cosmologically allowed lighter stop masses are shown. Large stop masses 
at moderate gluino masses are allowed due to Z-poles and Higgs-poles. The dark matter 
condition sets an upper bound of 1200 GeV for the lighter stop at large gluino masses. In 
the lower left picture the lighter stop masses are shown for the NMSSM. In contrast to the 
MSSM |[43|, the lighter stop mass is nearly proportional to the gluino mass as discussed 



in 



20| following from the suppression of the domain tuq \mi/ 2 \ in the NMSSM. In 
addition the constraints for the trilinear couplings A Ui , A^. and A e . (Eqs.(p^)- ([19|) ) are 
stronger in the NMSSM than in the MSSM. The corresponding constraints in the MSSM 
are obtained by replacing m 2 H in (Eqs.(|IT|)-(|T9"|)) by m 2 H +/i 2 . The stronger constraints 
for the trilinear couplings lead to a larger mass of the stop in the NMSSM, since the stop 
mixing depends on A t . Together with the UFB-3 constraint (Eq. fl27D ), which excludes 
small values of m , the lower bound for the lighter stop increases to 240 GeV in the 
NMSSM. The lower right picture shows the cosmologically restricted solutions. The lower 
bound for the lighter stop slightly improves to 270 GeV. The cosmological upper bound 
of 1300 GeV for the lighter stop in the NMSSM is similar to the corresponding upper 
bound for the lighter stop in the MSSM at large gluino masses. 

In the upper left picture of Fig.|8| we show the mass of the lighter scalar Higgs Si versus 
| tan j3\ for the MSSM. In the MSSM the mass of the scalar Higgs S\ has to be larger than 



the lower experimental bound of 64.5 GeV p0 |. The upper bound for ms x is 155 GeV . 



For | tan /3| there is a lower bound of about 1.4. In the upper right picture the dark 
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matter condition is imposed improving the upper bound for ms x to about 145 GeV . The 
lower left picture shows the mass of the lightest scalar Higgs Si in the NMSSM, which is 
lower than 150 GeV. Before applying the UFB-3 constraint (Eq. (|27j) ) there exist scalar 
Higgses Si in the NMSSM lighter than the MSSM Higgs bound of 64.5 GeV [0, which 
are predominantly a Higgs singlet N || [II], [E|. These very light scalar Higgses are now 
excluded by the UFB-3 constraint (Eq.(^7])), which excludes the parameter region with 
small m and large |mi/ 2 |. But this parameter region is needed in order to have very light 
scalar Higgses with a large Higgs singlet portion fT3|. In contrast to the MSSM, the lower 



bound for | tan/3| in the NMSSM is about 2.4. The higher lower bound for | tan/3| follows 
from the extra third minimum condition (Eq. ([TT|)) in the NMSSM. For larger values of 
| tan /3| the density of NMSSM solutions decreases, if the UFB-3 constraint (Eq.(|27|)) is 
applied. This follows from the fact that the NMSSM solutions with larger |tan/3| are 
lying in the parameter region with small mo and large |mi/2| H lOH - The lower right 



picture shows only solutions, which fulfill the dark matter condition. In the NMSSM the 
cosmologically allowed scalar Higgs Si is heavier than 69 GeV and lighter than 130 GeV. 
In contrast to the MSSM the dark matter condition together with the UFB-3 constraint 
(Eq. (|27|) ) leads to an upper bound of about 9 for |tan/3| in the NMSSM. However our 
investigations are restricted to | tan/3| values below 20 as mentioned in section 3. 

In the upper left picture of Fig.|9| we show the mass of the lighter scalar Higgs Si versus the 
mass of the pseudoscalar Higgs Pi for the MSSM. In the MSSM mp 1 is larger than 150GeV 
within our restricted range of |tan/3| values below 20 as mentioned in section 3. The 
upper right picture shows the cosmologically allowed Higgses. The dark matter condition 
restricts the mass of the pseudoscalar Higgs Pi to the range 200GeV < mp 1 < 2300Gel^. 
The lower left picture shows the mass of the lightest scalar Higgs Si versus the mass of 
the lighter pseudoscalar Higgs Pi for the NMSSM. If the UFB-3 constraint (Eq.(|2"T|)) is 
not applied there exist very light scalar Higgses Si correlated with very light pseudoscalar 
Higgses Pi which have both a large Higgs singlet portion ||. With the UFB-3 constraint 
(Eq. (|27|)) these very light Higgses are now excluded as already discussed for the very light 
scalar Higgs Si (see Fig.[§p. For this reason we have now a lower bound of 200 GeV for 
the pseudoscalar Higgs Pi in the NMSSM for | tan/3| values below 20. In the lower right 
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picture the dark matter condition is imposed leading to a lower bound of 250 GeV and 
an upper bound of 900 GeV for the lighter pseudoscalar Higgs P\. Compared with the 
MSSM the upper bound for the lighter pseudoscalar Higgs P\ is somewhat lower in the 
NMSSM. This follows from the extra third minimum condition (Eq.(PD) in the NMSSM 
and the suppression of the domain mo ^> |mi/2| in the NMSSM. 



7 Conclusions 

In the present paper we have made a comparison between the MSSM and NMSSM con- 
cerning the parameter space and the particle spectra. In contrast to the MSSM, the case 
m ^> \mi/2\ is disfavored in the NMSSM as discussed in and the squark masses are 
nearly proportional to the gluino mass (cf. Figs. (6) and (7)). Furthermore the gluino and 
the lighter u-squark in the NMSSM are heavier than in the MSSM. In the NMSSM the 
lower bound for the gluino is 330 GeV and 370 GeV for the lighter u-squark, whereas the 
lower bound for the gluino and lighter u-squark in the MSSM is 250 GeV. In the NMSSM 
a light stop is excluded by the stronger constraints for the trilinear couplings A Ut , 
and A £i (Eqs.([T7|)-(|T9D) and the UFB-3 constraint (Eq.(pTD). This leads to a larger stop 
mass m% x > 240 GeV in the NMSSM even before considering dark matter constraints. 
This is in contrast to the MSSM where stops as light as the experimental lower bound of 
63 GeV are allowed. One important difference between the NMSSM and the MSSM 



is the occurrence of LSP singlinos and light scalar and pseudoscalar Higgs singlets in the 
NMSSM. Applying the UFB-3 constraint (Eq.fl2"7p) these singlet states are excluded as 
shown in Fig. 4 and Fig. 9 and discussed in section 6. In this way the UFB-3 constraint 
strongly reduces the differences between the NMSSM and the MSSM. 
Furthermore we have investigated the consequences of the dark matter condition for the 
MSSM and NMSSM. Together with the UFB-3 constraint (Eq.([27|)) the dark matter con- 
dition gives an upper bound of about 9 for |tan/3| in the NMSSM (see Fig. 8), whereas 
in the MSSM larger values of |tan/3| are furthermore allowed within our investigated 
range of | tan/3| values below 20. For the LSP neutralinos, lightest charginos, gluinos and 
scalar Higgses Si the cosmologically allowed mass ranges are comparable for the MSSM 
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(40 GeV < m k o < 240 GeV, 85.4 GeV < m-+ < 450 GeV, 250 GeV < m- g < 1600 GeV, 
64.5 GeV < m Sl < 145 GeV) and for the NMSSM (50 GeV < m^o < 260 GeV, 
90 GeV < m-+ < 500 GeV, 370 GeV < m g < 1800 GeV, 69 GeV < m Sl < 130 GeV). 
The enhanced lower bound of 50 GeV and the larger Higgsino portion of the lightest 
cosmologically allowed LSP neutralinos leads to improved lower bounds for the charginos 
and gluinos in the NMSSM. Furthermore the dark matter condition slightly improves the 
lower bounds for the lighter u-squark and lighter stop in the NMSSM to 390 GeV and 
270 GeV , respectively. In the MSSM, due to Z-pole and Higgs-pole enhanced neutralino 
pair annihilation, large masses of the selectrons, u-squarks and stops are cosmologically 
allowed. In contrast, the NMSSM does not allow for large slepton and squark masses at 
moderate gluino masses (cf. Figs. (5-7)) and the dark matter condition leads to the upper 
bounds rrie R < 450 GeV , m Sfl < 1600 GeV and < 1300 GeV at large gluino masses, 
which are comparable with the MSSM. In comparison with the MSSM the upper bound 
for the cosmologically allowed pseudoscalar Higgs Pi is somewhat lower in the NMSSM 
m Pl < 900 GeV (Fig.9). 

Despite of the discussed differences, it may be difficult to distinguish experimentally be- 
tween the MSSM and NMSSM. One possibility is the measurement of smaller couplings 
of the neutral Higgses in the NMSSM, which result from the mixing of the ordinary Higgs 



states with the Higgs singlets as discussed in || |§ and [Tl|] - |T3[- Another possibility 
suggested in || [T0| is the search for an additional cascade in the NMSSM, which is pro- 
duced in the decay of the bino-like second lightest neutralino into the LSP singlino. Since 
the LSP singlinos are excluded by the UFB-3 constraint this additional cascade is not a 
way to discriminate between the NMSSM and MSSM any more. 
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Appendix Al: Annihilation cross section 



To calculate the relic neutralino density (Eq.([J{])) in the NMSSM we need the annihilation 
cross section of the neutralinos into the different decay products X and Y. 

X%h) x\{h) -> X{X X ) Y(X Y ) 

Xf = Xx — Ay, where Ax and Ay are the helicities of particle X and Y. Within the partial 
wave formalism the annihilation cross section is obtained |34| 



as: 



Vann(Xl Xl "> X Y)v 



A8nsS 



\ACSo)\ 2 + l(\A( 3 P )\ 2 + |A( 3 P0r + |A( 3 P 2 )| 2 ) 



(34) 



g W!. 2 « + + (^-^) a (35) 
s s 



s pa 4to| is the center of mass energy squared. S is a symmetry factor, which is 2 if 



X — Y. Summation over final helicities is implicit in Eq . (p4|) . The partial wave ampli- 
tude A( 2S+1 Lj) describes annihilation from an initial state with spin S, orbital angular 
momentum L and total angular momentum J. The partial wave amplitudes A( 2S+1 Lj) 
for the NMSSM are obtained by slightly modifying the partial wave amplitudes of the 



MSSM in [34[] and by using the corresponding NMSSM couplings. Thereby we have cor- 



rected some misprints [44| in pl| |. One modification of the neutralino annihilation is the 



exchange of 5 neutralinos in the NMSSM instead of 4 neutralinos in the MSSM. Other 
modifications concern the exchange of the scalar Higgses Si, 5*2 and S3 in the NMSSM 
instead of Hi and H2 in the MSSM and the exchange of the pseudoscalar Higgses Pi and 
P 2 in the NMSSM instead of P in the MSSM. In addition the NMSSM coupling g a _ + 
(Eq.(|98|)) has to be inserted in Eq. (p6|) . In the MSSM there is only one pseudoscalar Higgs 
P and the coupling g a - + turns out to be 1. In the following the partial wave amplitudes 
A( 2S+1 Lj) of the NMSSM for all decay channels are given in leading order in the relative 
velocity v of the non-relativistic neutralinos. Previously we give some NMSSM couplings 
and definitions. The couplings of the W boson to the chargino xt an d the LSP neutralino, 
which is denoted with index in [|34| instead of 1, have the same form in the MSSM and 
NMSSM, but a different matrix N (Eq.fljjD) enters: 

= ~ N 0A V j2 + N 02 V 3 i (36) 
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O^ = ^=N Q3 U j2 + N 02 U j1 . (37) 
The mass matrix of the charginos is diagonalized by the orthogonal matrices U and V, 



which are explicitly given in [ 40 . In the NMSSM the couplings of the charged Higgs to 



the chargino xt an d the LSP neutralino has to be modified compared to the MSSM: 

Q% = N 04 V n + 4= (iV 02 + N 01 tan 6 W ) V j2 - iV 05 V j2 (38) 

V2 92 

Q'o* = Nos U jt -^= (N 02 + N 01 tan B w ) U j2 N 05 U j2 . (39) 

\/2 92 tanp 

In the MSSM and NMSSM the coupling of the Z boson to the neutralinos Xi an d X° is 
given by: 

0'[f = ~ N l3 N j3 + ~ iV l4 N H . (40) 

In the NMSSM the coupling of the scalar Higgs S a to the neutralinos x° and x°j changes 
compared to the MSSM and can be written as: 

Tsai^-U^Q'tj + UgS^ + U^Z^ (41) 
Q"j = 7T~ l N * (92 N j2 - g y N jX ) + y/2 X N l4 N j5 + (i ~ j)] (42) 

1 92 

S"j = 7^- [N* (92 N j2 - g y N n ) -V2\N l3 N j5 + (i <- j)) (43) 

^92 

Z'L = TT^hv^ A N t3 N j4 + V2 k N i5 N j5 + (i <-> j)). (44) 
2 g 2 

The 3x3 mass matrix of the scalar Higgses Si [[|] in the basis {H® R , H^, N R ) with 
Hf R = v2 Re{H®} (i = 1 — 3) is diagonalised by the matrix II s . The scalar Higgses 
Si (i = 1 — 3) are ordered with increasing mass. The coupling of the pseudoscalar Higgs 
P a to the neutralinos is also different in the MSSM and NMSSM. In the NMSSM this 
coupling has the following form: 

Tp ai j = —U^ Q"j + U^ 2 S"j -U^ 3 Z"- (45) 
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Qij = 7^- [Na (92 N j2 - g y N n ) - V2 A N l4 N j5 + (i <- j)} (46) 
Sij = 7f-Wn (92 N j2 - g y N n ) + V2XN t3 N j5 + (z <- j)] (47) 

^ 92 

Z'H = Z%. (48) 

The 3x3 mass matrix of the pseudoscalar Higgses P a and the neutral Goldstone boson 
§ in the basis (H^, H$ T , iV>) with = V2Im{H®} (i = 1 - 3) is here diagonalised by 
the 3x3 matrix £/ p . The pseudoscalar Higgses P a (z = 1, 2) are ordered with increasing 
mass. In ]34[] the following definitions are used: 



A 2 = m * + my (49) 
2m| V ; 

A. 



i? x = !M etc. (50) 

m x 



P I = l + Rj-(R 2 x + R 2 Y )/2. (51) 



i. xUi^w-(\)w + (\). 

The W bosons are produced by Z boson and scalar Higgs Si (i = 1,2,3) exchange in the 
s-channel and chargino x~j (j = 1> 2) exchange in the t- and u-channel. 



A^Sq) : X f = 0,A = ±1 



(OL2 _j_ f)R2 



+ V2v 2 /3 f g 2 2 



2 R] 



L nL qR , 

3 p2 ^Oj ^Oj + 



o L2 + o R2 



P; 



2 1 



3P j 



2 13 
3P 2 



(52) 



A( 3 P ):A/ = 0,A = 



v/6 



R2 



4 Off + Off 4QfcOgiff 
3 R R 



3P, 
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Otf + Oof 



3 Pi 



2 0feQg^ 



3R 



VEvjl + P^gtFi 
(A-m+iGsjRw 



0, A = ±1 



0^ + °6? 



1_2# 
3 3P, 



2 Og Og i?+ 



+ 



4 - ifc + iG Si 



Pi) ■■ |A/| = 1 



2vP f \ f gl 



R 



w 



off + o f 



1 - 



2 0£0g^ 
^ 2 



+ 2v/3 f gt 



QL2 _ QR2 
2 ^Oj ^Oj 



Ryy Pj 



2- 



8vP f g*Q>£ 
R w (4 — 



= 0, A = 



n W r 3 



(4 - R%) R 



= 0, A = ±1 



2v Pf n 2 (n L2 nR 2, ±v(3 f glO^ 

p #2 l^Oj ^Oj J 



4-i?| 



^):|A/|=2 



1vg\ 



R 



w 



~ J jf(Off + 0™) + 

3 



2 O^, O^- —2 



n L2 n R2 
Oj ^Oj lb j -75 2 , x 75 3 U 0i ^Oj 
2 />/ + A / />/ 52 

/' 3 
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A/ = 0, A = ±1 



^_ ^cy 0j -t- cy 0i j 



1 -P 



+2 



R 2 



(60) 



\ f = o, A = 



4^1 
v^P 2 



^ 2 + O f 



9 7? 2 1 



+ 4/3' 



-2 Oqj O 0j '"'J 



^ 2 



(61) 



G*5 i = r s . m S - / m| are the rescaled widths of the scalar Higgses The coupling of the 
scalar Higgs to the LSP neutralinos and to the W bosons is combined to Pj. In the 
NMSSM F t is given by: 



F i = (Uf 1 cos + Uf 2 smP)T Si00 . 



(62) 



2. x?xS|-Z(A)Z(A). 

The Z bosons are produced by scalar Higgs Si (i — 1,2, 3) exchange in the s-channel and 
neutralino x° (j — 1 — 5) exchange in the t- and u-channel. 



A^Sb) : A/ = 0, A = ±1 

4v / 2^ / ^ 2 O / f 2y/2v 2 P f g%0. 



cos 2 # w -Pj 
A( 3 Po):A / = 0,A = 



+ 



//L2 



cos 2 9 W 



cos 2 # w R 2 7 0j 



Rj 1 



4_ 2fi 
" + 3P 2 



3 P. 



«L2 



Pi 



3P P 



3P 



(63) 



+ 



cos 2 0^ 



1 - 



3P, 



+ 



Pj 



1 - 



3P 



(4 _ tf^ + iGsi) Rw 



(64) 



A/ = 0, A = ±1 



2y/6^ 2 O f 
cos 2 ^ w 



1 




P, 


5 


3 3P, 


H 

P 

r 3 



VEvg 2 . R w Ft 



cos 2 ^ (4-R 2 Si +iG Si ) 



(65) 
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M 3 Pi) : |A/| = 1 



4vP f \ f giO% 
cos 2 9 W R z 



I/L2 



i-l 



Ri 

-I i. 



(66) 



^2) : |A/| = 2 



4^ 



COS* v w 



&w Pi 



pill LI 
W 0j 



(67) 



|A/| = 1 



G „L2 

cos 2 e w R z 0j 



+ 



3 



f 



(6S 



\ f = 0,A = ±1 



Avgj 



cos 2 9w 



r\liL2 
U 0j 



E>2 E>2 

- Kj - itg 



(69) 



A/ = 0, A = 



^3 COS 6yy Ry 



f)llL2 
U 0j 



R%P 



■2 1 



+ 



2 Rip, 



(70) 



3. xlxl^Z(X)S a . 

The Z boson and the scalar Higgs S a are produced by Z boson and pseudoscalar Higgs 
P a (a — 1, 2) exchange in the s-channel and neutralino x° (j — 1 ~~ 5) exchange in the t- 
and u-channel. 



A^So) : A = 



2y/2P f 



cos 9 



w 



2J 3 {R 3 -l) 



+ 



O 



Pj ' Rz cos 9w 4 



4 Lg 

R 2 Pa + iG Pa 



— v 



V2P f 



9l 



J. 



Ri 



cosdw 



(Rj - 5) - 



2(^-1) , 4(^-1) 
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3P 2 



P f + (2-A 2 
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+ V 



Rz cos 9} 



3L r 



O 



A-Rp+iG Pa 4 R z cos 6 W 



(71) 



A^P,) : A = ±1 



9l 



J; 



4 V ^A L 

COS^^y P. 2 



8 A 
3 4 



g\ Jj Rj 



cos 9 W P. 



+ 2v 



92 q R z 



cos 2 9 W " 4 - P 2 



(72) 



A = 



-2 



P; 



i + 5 - 

2 



S3 



COS6 1 



2(l + ^)#- 



P; 



P, cos^(4-P|) 



(73) 



A( 3 P 2 ) : A = ±1 



1 costV Pf 



In [B4| the following definition is used: 



5' = (R z -Ri)/2. 



(74) 



(75) 



Gp a = Tp a mp a /m^ are the rescaled widths of the pseudoscalar Higgses P a . The couplings 



appearing in the neutralino Xj exchange are combined to Jf 

Jj = -0'of T SaOj- 



(76) 



L a is the product of the coupling of the pseudoscalar Higgs P a to the LSP neutralinos 
and to the Z boson and the scalar Higgs S a . In the NMSSM L a reads: 

1 



La — ^ ( U a2 Ua2 ~ ^al ^al) Tp a 0O- 



(77) 



O consists of the couplings appearing in the Z boson exchange. In the NMSSM O is given 
by: 



= 0^([/f lC os/3 + [/* sin/3). 



(7S 
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4. x? jg -> Z(A) P a . 



The Z boson and the pseudoscalar Higgs P a are produced by scalar Higgs Si (i = 1,2,3) 
exchange in the s-channel and neutralino x'j ( J = 1 ~ 5) exchange in the t- and u-channel. 



A( 3 P ) : A = 



Ay/6 



Rz cos 6 



2 



3Pf 



2 - 2 A 2 S 21 
- + R 2 h - 

3 J 3 6 



*k 4 L » 
P 2 A 

3 



R Si 



(79) 



A( 3 Pi) : A = ±1 



AvP f X 



9l 



COS 9w 



-R3 + 



*0) 



A( 3 P 2 ) : A = ±1 



1 cos 6 



w 



s21 



-R j + 



P 2 

j 



(81) 



A = 



v(3 



15 92 



^3 1 cos6 w R Z P 2 



1 + Rj - A 2 + — (Rj - 1) 



(82) 



In [p4j the following definition is used: 



6 2 = (Rl - R P J/2. (83) 

The couplings appearing in the neutralino Xj exchange are combined to Kf 

K s = OgT Pa0j . (84) 

L{ is the product of the coupling of the scalar Higgs Si to the LSP neutralinos and to the 
Z boson and the pseudoscalar Higgs P a . In the NMSSM Li has the following form: 

1 



L% = 2 ( Ui2 U » 2 ~ U K U ai) TsiOo- 



5) 
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5. xlxl^W-{\)H + . 



The W boson and charged Higgs are produced by scalar Higgs Si (i — 1, 2, 3) and pseu- 
doscalar Higgs P a (a = 1,2) exchange in the s-channel and chargino Xj (j — 1) 2) exchange 
in the t- and u-channel. 



ACS ) : A = 



r - „J' + J" r 2 f3 f g 2 2 (J' + J") 

rt\Y rj 



Rw Pj 



, -7$ 



5 2 4/^ 2 2 
2"^ + 3^"3^ ( ) 



3 J? P,- + 2 



- AV2p f g 2 2 g a _ + T Pc 



1 



oo 



i?w 4 



1 + 



(86) 



A( 3 P ) : A = 



*7 



>4 



3P/ 



- + i?f --A 2 + - 

3 J 3 6 



*2 



4v/6 -3 



R 



w 



Si 



(87) 



A( 3 Pi) : A = ±1 



— Av 



S 4 

J 4 



7' 4- J" R + 

r 3 r 



AvdrX 9 

+ — ST - ^ 

i 



Rl{K>-K>>)--(Ki + K»\ 



(88) 



A = 



>4 



(89) 
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A( 3 P 2 ) : A = ±1 



^ 2 



^ 2 



^2 



(90) 



A = 



V3R w Pf 



C2l 



-1 + A 2 + 



1 + 



*7) 



(91) 



In [[34] the following definition is used: 



6 2 = (R 2 w -R 2 H+ )/2. 



(92) 



The couplings appearing in the chargino xt exchange are combined to: 



■3 



4 (0§ - Oo;) « sin/? + Qg cos/5) 



(93) 



= "J (0§ + (Og sin/3 - Qg cos/5) 



(94) 



*5 



(0§-0§) (Qgsin/3-ggcos/3) 



(95) 



*7 



J (Og + 0§) (Qg sin/5 + gg cos/3). 



(96) 



Lj is the product of the coupling of the scalar Higgs Si to the LSP neutralinos and to the 
W boson and the charged Higgs. In the NMSSM Li is given by: 

1 



Li = ^ (^2 cos/5 - tf£ sin/?) T 5j00 . 



(97) 



In the NMSSM the coupling g a (_ of the pseudoscalar Higgs P a to the W boson and 

charged Higgs reads: 



g a - + = 1/^2 cos (3 + U*i sin/3. 



(9* 



In the MSSM with only one pseudoscalar Higgs P the coupling g a h turns out to be 1 |34 
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6. Xi Xi - Sa S b or P a Pp. 



The scalar Higgses S a and S b are produced by scalar Higgs Si(i = 1,2,3) exchange in the 
s-channel and neutralino Xj ( J = 1 — 5) exchange in the t- and u-channel. 



A( 3 P ) : 



9abi T t 



R> 



Si 00 



4 - R 2 + iG Sl 



-~2 ■ 

OT T 1 + R 3 , 4 T T 

* ± SaOj ± S b 0j — p r — J SaOj J- S b 0j p2 



(99) 



A( 3 P 2 ) 



~ v Pf 9 2 Ts a oj T Sb0j pa 



(100) 



The pseudoscalar Higgses P Q and Pp are also produced by scalar Higgs Si (i = 1,2,3) 
exchange and neutralino y*j ( J = 1 — 5) exchange. The corresponding amplitudes for P a 
and Pp production can be obtained by replacing a, b by a, f3 or S a , S b by P a , Pp and Rj 
by — Rj . In the NMSSM the trilinear scalar Higgs coupling g aaa is given by: 

6 



9a 



92 ra z 



9y + 9l 



( Vl U S J + v 2 U S J - Vl U s al Ug 2 - v 2 Ug 2 Ug 2 ) 



+^= 2 V2 ug 2 ug 2 + x ug 3 ug 2 + ^= Vl ug ug 2 + -| . ug 3 

+^ «i ^ 2 + ^ ^ ^ + v/2 A; 2 x 
- «i ^ u!i ~^v 2 Ug Ug 2 -yftXkxUZ U s a2 Ug 3 



A A 



A r/ 5 u s u s - 

U al U a2 U a3 



kAi 



U a3 



V2 ~ a ^ a ^°* 3^2 
The trilinear scalar Higgs coupling g aab (or g aba ) has the following form: 



(101) 



9aab 9aba 



92 m z 



9 -^-^ V2 (3 Vl Ug 2 + 3 v 2 Ug Ug 2 - Vl Ug Ug 2 



-2 Vl Ug Ug U s a2 - v 2 Ug Ug, 2 - 2 v 2 Ug 2 Ug Ug,) 
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^X 2 v 2 - 



V2 XkVl (TjS TjS2 , 2 tjS tjS tjS \ , /rrS rrS2 , 9 rrS rrS rrS \ 

^2 v 62 a3 + a2 63 a3 / + ^2 ^ 3 2 a3 62 a2 > 

X 2 v 1 -Xkv 2 . s S2 s s s , A 2 ^ . s S2 s s s , 

^2 \ U bl U a3 + ZU al U b3 U a3) + ^ \ U b3 U al + ZU a3 U bl U al) 

^1 (rrS TjS2, 2 rrS rrS rrS \ , ^ ^ 2 (rrS rrS2 , 9 rrS rrS rrS \ 



-f 

A 



6fc^rfc4 ..5 ,,g 2 

~^ ^7f 3 3 



2 A kx + A A 



A K U s a2 Ug + t£ tfg C/f 3 + Ug 2 Ug)^j . (102) 



v/2 

Finally the trilinear scalar Higgs coupling gi 23 (or #i 32 or g 23i ) reads: 



S123 = S132 =9231 = — ( V2(6 Vl Ug Ug Ug + 6 v 2 Ug C/| 2 f/f 2 

#2 V ° 

-2 Vl (ug u 2 s 2 u! 2 + Ug ug u! 2 + ug ug U 2 S 2 ) 

-2 v 2 (Ug [/£ [/£ + [/£ E# Ug + 0* [/£) ) 
+ ( A 2 v 2 - A fc Wl ) (C/f 2 Ug Ug + tf£ C/f 3 Ug + C/ 3 5 2 C/f 3 C/| 3 ) 

+V2X 2 x (Ug u 2 s 2 ui 2 + ug ug u! 2 + ug ug U 2 S 2 ) 

+V2(X 2 v 1 -Xkv 2 ) (Ug Ug Ug + Ug Ug Ug + Ug Ug Ug) 

+V2X 2 x (ug ug ug + ug ug ug + ug ug ug) 
+V2X 2 Vl (ug ug ug + ug ug ug + ug ug ug) 
+V2X 2 v 2 (ug ug ug + ug ug ug + ug ug ug) 

+V2 (Qk 2 x-k A k ) Ug Ug Ug 

( II s II s II s 4- II s II s II s 4- II s II s u' 

\ u ll u 22 u 33 T <^ii ^32 u 23 ^ U 21 u 12 u 3 



2 Xkx + XA X , s s s s s s s s s 

^2 V ! l ( 22 ( !3 ' L J. 1 32 ' 23 ' 2 : 1 J 2 ' 33 



-ug ug ug + [/£ ug ug + [/£ f/ 2 5 2 ug) ^ . (103) 



The trilinear scalar and pseudoscalar Higgs coupling g aaa (or g aa a) is given by: 

,2 i ^2 



gaaa = 9aaa = — ( ^ («1 ^1 ^ + ^2 C/f 2 U t 

92 mz V ° 



P2 
a2 
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-vi u s al u p J - v 2 U s a2 U P2 ) + ^ (v 2 U s a2 U P J + x u P2 + Vl u s al U P J 

U P J + v, U s al U P2 + v 2 U s a2 U P2 ) 
+ V2 k 2 x U s a3 U P J + ^ ( Vl U s a2 U P J + „ 2 t£ U P J) 
-V2\k{ Vl U s a3 C/ p t£ + V2 [/f 3 t£ E£ 
t£ t£ t& + x U s a2 U p aX U p 3 - x U s a3 U p al U p 2 ) 

+ ^ m u p 2 u p 3 + u s a2 u p al u p 3 + e£ u p al u p 2 ) + ^ e£ • (104) 

The trilinear scalar and pseudoscalar Higgs coupling g a p a (or g aa p or g a /3a) can be written 



as: 



1 ( 9y 92 / — / P P 

9af3a = 9aaf3 = 9a(3a = " V 2 ( Vl U al U al U 01 

92 f^z \ o 
+v 2 U s a2 U p 2 U p 2 - Vl U s al U p 2 U p 2 - v 2 U s a2 U p al U p x ) 

+ y/2 X 2 {v 2 U s a2 U p 3 U p 3 + x U s a3 U p 2 U p 2 + Vl U s al U p 3 U p 3 

+x U s a3 U p al U p x + Vl U s al U p 2 U p 2 + v 2 U s a2 U p aX U p x ) +2V2k 2 x U s a3 U p 3 U p 3 

+V2\kv 1 U s a2 U p 3 U p 3 -V2\k Vl U s a3 (U p 3 U p 2 + U p 3 U p 2 ) 

+V2Xkv 2 U s al U p 3 U p 3 -V2Xkv 2 U s a3 {U p 3 U p x + U p 3 U P A 
-2 A k x + A A x ( s . P P P P , s , P P P P 

U al \ U a2 U /33 + U (52 U a.3) + U a2 \ U al ° 133 + U f31 U a3 



+ 2Xkx J_ XA " U s a3 {U p al U p 2 + U p x U p 2 ) +V2kA k U! 3 U p 3 C/ p ) • (105) 



7. xlxl^S a P p . 



The scalar Higgs S a and the pseudoscalar Higgs Pp are produced by Z boson and pseu- 
doscalar Higgs P a (a = 1,2) exchange in the s-channel and neutralino x°j (j ' — 1 — 5) 
exchange in the t- and u-channel. 



AfSo) : 



2V2g 2 2 



Rz 



A-R 2 P 



+ V2g 2 2 9a/3Z 



K P ~ K S a 

R\ 



l- v - 

8 
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R ; 



-4V2g 2 2 T Sa0j T Pp0j -L 



l + v 2 



8 2R 



i + Ft 



w 2 



\/2 (R 2 P - R 2 S J g\ T Sa0 j T P ^ 0j 



1 + v 2 



8 IP. 



i + ft 



3P 2 



(106) 



A( 3 Px) : 



92 9a/3Z 



U-R 2 + iG 2 



Av(3 f g 



(107) 



Gz = T z m z I m| is the rescaled widths of the Z boson. The couplings appearing in the 
Z boson exchange are combined to g a pz- In the NMSSM g a/3Z has the following form: 

1 



9a/3Z — 



2 cos 2 9w 



m u; 2 - U s al U p pi ) 0#. 



(108) 



8. x\x"^H + H-. 

The charged Higgses are produced by Z boson and scalar Higgs Si (i = 1,2,3) exchange 
in the s-channel and chargino Xj C? = 1> 2) exchange in the t- and u-channel. 



A( 3 P ) 



-\fcvgl 



(Q'of + Q 0j 



1 2fi 



3P 2 



R 

2 Qoj Qoj 



A( 3 P 1 



A( 3 P 2 ) 



+ V6 v gl R w ^ + * 



4-^4 



(109) 



2 v ]3 f g\ (Q>? - Qg>) -L + 4 „ 0, r2 



W) „ 2 



cos 2 ^ ^4-i?| 



(110) 



-^vp 2 f9 UQ'o? + Q'o?)± 



(111) 
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g + -i is the product of the coupling of the scalar Higgs Si to the LSP neutralinos and to 
the charged Higgses. In the NMSSM g + _i is given by: 

cos(2/3) 



9+- 



Si 00 



(E/£sin/3 + E/£cos/3) + 



2 cos 2 9w 



-U* cosP + U? 2 sin (3) 



V2X 2 v 
92 m w 



(Uf 2 sin (3 cos 2 /3 + U^ sin 2 (3 cos (3) + 



^2A< 



92™W 



H E/£ sm /3 cos /3 H E/£ sin /3 cos /3 

92 mw 92m w 



(112) 



9. X? *?->/.(*) /«(*)■ 



The fermions are produced by Z boson, scalar Higgs Si (i = 1,2,3) and pseudoscalar 
Higgs P a (a = 1, 2) exchange in the s-channel and sfermion / 12 exchange in the t- and 
u-channel. 

A( l S ):X f = 



V2(-ir^(X'l + W' 2 a0 ) 



1 + v' 



1 + # 



2 Pi 3P! J 



5t 
p 



+2v / 2(-i) I+1 / 2 x; w/; -L 



1 + *r 



4 2 Pi 6 Pi 3P{ 



HKo - ^o, <o - Ko, Pi - P 2 ) + (-lf +1/2 2 ~^f O" L 00 T 3a |> 

COS u\\r tt z 



+ AV2(-lf +1 / 2 g 2 h Paa T Pa00 



4 - R 2 Pa + *G Pa 



1 + T 



(113) 



A( 3 P ) : A/ = 



1 

P 



3P 2 



+ fv]3 f (X>l + W>l )^ 



+ (X' a0 <-> z' a0 , w' M <-> r a ' , p <-> P 2 ) - 2 V6 v /?/ 0, 



PgjQO 

4 - P| + iG* 



(114) 
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A^P,) : \ f = 



^(X' 2 a0 - W ,2 a0 ) + (X' a0 <-> Z^o, <-> F a ' 0) P 1 <-> P 2 ) 



-2v9l^-[T 3a -2e fa sm>e w ] H " 



cos 2 9 W 



4-RI + iG z 



X f = ±1 



(115) 



V2v 



X f P f (X' 2 a0 + W' 2 aQ ) 



1 1 

~1\ + pI 



+ {X'lo - W' 2 a0 ) 



Pi pi 



+2(3 f X f X' a0 W^^ 



+ (X' o »"K ,K "Y^P 1 ~P 2 ) 



+ 2V2v^—O"% [X f T 3a P f -(T 3a -2e fa sin 2 9 W )] 



4-RI + iG z 



(116) 



A( 3 P 2 ) : A, = 



2 - 1 



R l {X'l B + W'l 1 )+2X' M W'„ 



aO 



+ (x' a0 <-> z' a0 , w' a0 <-> r a ' 0) Pi <-> p 2 ) 



(117) 



X f = ±1 



1 



« 0/ ^ [-(^'ao + W /2 a0 ) + (3 f Xf (X' 2 a0 - W' 2 a0 ) - 2 X^o W'J 



+ (x' »"Ko,K "Y^p 1 ~p 2 ). 



(118) 



The sfermion mass eigenstates /i )2 are defined by 



/i = /l cos 6^ + f R sm6 f 



(119) 



/2 = -h eon Of + / fl cosOf. 



(120) 



The couplings of the mixed sfermion /i )2 to the fermion and LSP neutralino can be written 
as: 



X' a0 = X a0 cos9 f + Z a0 sin 9 f 



(121) 
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W a0 = Z a0 cos9 f + Y a0 sin0/ 



(122) 



Z' a0 = -X a0 sin9 f + Z a0 cos6 f (123) 

Y; o = -Z a0 sm9 f + Y a0 cos9 f , (124) 

where the couplings of the unmixed sfermion to the fermion and LSP neutralino are 
given by: 

X a0 = -V2g 2 [T 3a N 02 - tan 9 W (T 3a - e fa ) N 01 ] (125) 
Y a0 = V2g 2 tan%e /a 7V i (126) 



^o = - / -' 2m ; iV 04 (127) 
V2 sin p 



Z« = - N 03 . (128) 

V 2 cos p mw 

In the NMSSM the couplings of the Higgses to the fermions have to be modified compared 
to the MSSM: 

g 2 m u Ug g 2 m d Ug 

2 sm p m w 2 cos p m w 

, 92 m u U? 2 g 2 m d 

2 sm p m w 2 cos p m w 
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Appendix A2: Renormalisation Group Equations 



For completeness we also give here the 1 loop SUSY RGEs for the parameters in 
the NMSSM with our conventions. Q is the running mass scale in the variable t = 
\og(Q/M weak ). 

Tt« = (131) 

dir. , 

dt* = m^ 91 (132) 

Tt 9 * = (133) 

= ( k 2 + 2A 2 + -h 2 - -g\ - -g 2 

dt 8vr 2 V 2 * 2 y2 2 yy 

8 7T Z 



m 2 



dt 

4t a * = ^(4\ 2 A x -2k 2 A k -3h 2 t A Us 

at O 7T Z 

- 3(? 2 M 2 - gjMx) (137) 
^ = (-\ 2 A X + k 2 A k ) (138) 

= (- A 2 A A + 3(1 + fe)^4 3 

+ y (?3 2 M 3 + 3s 2 M 2 + y^M^ (139) 

-r- ^r>- = — ( - X 2 A x + 5 i3 hf A v . 
dt 8vr 2 V 

+ y<7 3 2 M 3 + 3£ 2 M 2 + -gjMj (140) 
A* = ("A 2 A A + 3^ 2 M 2 + 3^Mx) (141) 

rfn 1 



dt ^ " 8tT 2 



[2 A 2 ( m 2 N + + m 2 ^ + A 2 ) 



+ 2k 2 (3m 2 N +A 2 k )] (142) 
1 

(//-'"'" 8~^ 



3g 2 2 M 2 - £ 2 M 2 ] (143) 

2 J_ 



m k = T7 m k + ("4 + m 2 U3 + m 2 Qa + A 2 Ua ) (144) 
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8tt 2 



<*« h 2 t ( m^ 2 + m 2 U:j + m 2 Q:i + A 2 Ua ) 

1 



8tt 2 



y glMt - 3^ 2 2 M 2 2 - ^M 2 



2 5 l3 fc? ( + m 2 u 3 + ™ 2 Qs + ^b, ) 
16 , , 16 



- — ^ Mf 



8tt 2 
1 

8tt : 
1 

8^ 
11 

8^ 

1 

8^2 



16 

y 



9l Ml 



2 (-Sg 2 2 M 2 - g 2 M 2 ) 



-4g 2 y M 2 ) 

g 2 M 2 



8tt 



2 g\ m 3 
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Figure Captions 



Fig 1 The scalar mass m versus the gaugino mass |toi/2|. The upper pictures are for 
the MSSM, the lower pictures are for the NMSSM. In the right pictures the dark 
matter condition is imposed. 

Fig 2 The trilinear coupling \Aq\ versus the scalar mass itlq. The upper pictures are for 
the MSSM, the lower pictures are for the NMSSM. In the right pictures the dark 
matter condition is imposed. 

Fig 3 The gaugino mass \M2\ versus the Higgs mixing parameter /i both taken at the 
electroweak scale. The upper pictures are for the MSSM, the lower pictures are for 
the NMSSM. In the right pictures the dark matter condition is imposed. 

Fig 4 The bino portion of the LSP neutralino versus the mass of the LSP neutralino. 
The upper pictures are for the MSSM, the lower pictures are for the NMSSM. In 
the right pictures the dark matter condition is imposed. 

Fig 5 The mass of the lighter chargino versus the lighter selectron mass. The upper 
pictures are for the MSSM, the lower pictures are for the NMSSM. In the right 
pictures the dark matter condition is imposed. 

Fig 6 The mass of the lighter u-squark versus the gluino mass. The upper pictures are 
for the MSSM, the lower pictures are for the NMSSM. In the right pictures the dark 
matter condition is imposed. 

Fig 7 The mass of the lighter stop versus the gluino mass. The upper pictures are for 
the MSSM, the lower pictures are for the NMSSM. In the right pictures the dark 
matter condition is imposed. 

Fig 8 The mass of the lightest scalar Higgs Si versus | tan/3|. The upper pictures are for 
the MSSM, the lower pictures are for the NMSSM. In the right pictures the dark 
matter condition is imposed. 
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Fig 9 The mass of the lightest scalar Higgs Si versus the mass of the lightest pseudoscalar 
Higgs Pi. The upper pictures are for the MSSM, the lower pictures are for the 
NMSSM. In the right pictures the dark matter condition is imposed. 
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Fig. 6 
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Fig. 9 
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